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c^ ■ We analyze the finite temperature behavior of the Sakai-Sugimoto model, which is a 

holographic dual of a theory which spontaneously breaks a U{Nf)L x U{Nf)ji chiral fiavor 
symmetry at zero temperature. The theory involved is a 4-1-1 dimensional supersymmetric 
SU{Nc) gauge theory compactified on a circle of radius R with anti-periodic boundary 
conditions for fermions, coupled to Nf left-handed quarks and Nf right-handed quarks 
which are localized at different points on the compact circle (separated by a distance 
L). In the supergravity limit which we analyze (corresponding in particular to the large 
Nc limit of the gauge theory), the theory undergoes a deconfinement phase transition 
at a temperature T^ = l/27cR. For quark separations obeying L > Lc — 0.97 * R the 
chiral symmetry is restored at this temperature, but for L < L^ — 0.97 * R there is 
an intermediate phase which is deconfined with broken chiral symmetry, and the chiral 
symmetry is restored at T^sb — 0.154/L. All of these phase transitions are of first order. 
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1 Introduction 

QCD, the theory of the strong nuclear interactions, exhibits two interesting phenomena 
at low energies - confinement and (approximate) spontaneous chiral symmetry breaking. 
Both of these phenomena are strong coupling effects which are not visible in perturbation 
theory, and there are no known quantitative methods to study them (except by lattice 
simulations). Apriori there is no relation between these two phenomena, and there are 
known examples of chiral symmetry breaking in non-confining theories and of confinement 
without chiral symmetry breaking. However, in QCD both phenomena occur together. 

New methods for studying strongly coupled gauge theories have been developed in the 
last decade, following the AdS/CFT correspondence [1]. These methods are particularly 
powerful in large Nc gauge theories which are dual to weakly coupled string theories in 
weakly curved spaces, since then many computations can be explicitly performed. This 
class of theories does not include (large Nc) QCD, but there are several known cases of 
theories which are labeled by a dimensionless parameter, such that in one limit of the 
parameter the string background is weakly curved, while the other limit gives (large A^^) 
QCD [21 El in El- It is believed that these theories are (for any value of the dimensionless 
parameter) in the same universality class as QCD, so it is interesting to try to use 
these theories to understand better the phenomena of confinement and chiral symmetry 
breaking. 

Understanding confinement in the weakly curved limit is amazingly easy; the confining 
strings are explicitly seen as strings in the weakly curved background (either fundamental 
strings or strings arising from D-branes), and their tensions are easily computable. It 
would be nice to be able to similarly study the phenomenon of chiral symmetry break- 
ing (xSB), and to see if (in the weakly curved limit) there is any link between x^B 
and confinement. Unfortunately, it seems much more difficult to find dual models for 
theories with non-Abelian U{Nf)L x U{Nf)fi chiral symmetry than to find dual models 
for confining theories. In the context of holographic dual models, quarks in the fun- 
damental representation of SU{Nc) were first incorporated in terms of probe D branes 
in jH]^. Since then several holographic models dual to theories with quarks have been 
constructed and investigated, with the quarks coming from additional "fiavor D-branes" , 
in both conformal and confining gauge theories [HI O UHl UH 1121 UHl UH UHl- However, 
there are relatively few brane constructions with non-Abelian chiral symmetry (see, e.g., 
|16j). and it has proven difficult to use these brane constructions to construct (by taking 
a near-horizon limit) string theory duals of theories with non-Abelian chiral symmetry 
(with or without xSB)^. 

^The first appearance of quarks in the fundamental representation of Sp{Nc) was in [?]• 

^On the other hand, there are several models realizing a C/(l) axial flavor symmetry, such as [nillOj. 



This situation was recently changed by the work of Sakai and Sugimoto fZlE], who 
found a simple brane construction realizing chiral symmetry, and took its near-horizon 
limit to find a string theory dual to a field theory with chiral symmetry breaking. The 
brane construction is in type IIA string theory compactified on a circle of radius R 
in the X4, direction, and involves Nc D4-branes wrapped on the circle (filling the 01234 
directions), Nf D8-branes sitting at the point 0:4 = and Nf anti-D8-branes sitting at the 
point X4 = L ^. The circle could have either periodic or anti-periodic boundary conditions 
for the fermions; the more interesting case for relating to QCD is the case in which all 
fermions have anti-periodic boundary conditions around the circle. The gauge theory on 
the D4-branes at low energies is a U{Nc) 4 + 1 dimensional maximally supersymmetric 
Yang-Mills theory, which is coupled to Nf left-handed 3 + 1-dimensional chiral fermions 
(in the fundamental representation of U{Nc)) at the intersection with the D8-branes and 
to Nf right-handed chiral fermions at the intersection with the anti-D8-branes. With 
anti-periodic boundary conditions, the fermions in the 4 + 1 dimensional gauge multiplet 
obtain a mass of order 1/R, and the scalars obtain masses from quantum corrections, 
so the (classical) massless spectrum is that of a 3 + 1 dimensional U{Nc) gauge theory 
coupled to Nf massless Dirac fermions. The U{Nf)L x U{Nf)ii global chiral symmetry 
of this theory^ is explicitly visible as the gauge symmetry on the D8 and anti-D8-branes. 

The string theory described above is non-supersymmetric and unstable, but this insta- 
bility does not affect the low-energy behavior, and disappears when we take a decoupling 
limit for the theory on the D4-branes; the unstable modes, which are the scalar fields cor- 
responding to the radius of the circle and to the asymptotic brane-anti-brane separation, 
become non-normalizable in this near-horizon limit so we do not need to worry about 
them. This decoupling limit, before the compactification on a circle, was described in 
[T^ , and it involves taking the string scale a' to zero and the asymptotic string coupling 
gs to infinity while keeping fixed the 4 + 1 dimensional gauge coupling g^ ~ gs^/a'. The 
fact that the (asymptotic) string coupling goes to infinity in this limit is related to the 
fact that the 4 + 1 dimensional gauge theory is not well-defined at high energies, and 
should be defined with some UV completion; in type IIA string theory the completion 
defines this theory as the 5 + 1-dimensional A/" = (2, 0) superconformal theory of M5- 
branes compactified on a circle, but this will not be relevant for our purposes (since we 
will always be interested in energy scales well below l/f^f). The near-horizon limit of 
the D4-branes, which is dual to their (UV-completed) 4 + 1-dimensional supersymmet- 
ric gauge theory, was described in [19^. It is highly curved at small radial coordinates 
(corresponding to the fact that the gauge theory is IR-free), and at very large radial 

^In 22] the specific choice L = ttR was made, for which the D8's and anti-D8's are antipodal on the 
circle, but this is not really necessary. 

^The axial U{1) is anomalous, but this is a subleading effect in the large Nc limit which will not be 
relevant for our leading order computations. 



coordinates the string coupling becomes large and the background needs to be lifted to 
M theory, but for large A*";, there is a large range of values of the radial coordinate for 
which type IIA string theory gives a good description. 

The near-horizon limit of D4-branes compactified on a circle with anti-periodic bound- 
ary conditions was first discussed in [2], where it was noted that it was the same as the 
(Wick rotation of the) near-horizon limit of near-extremal D4-branes, and that when the 
radius of the circle is small compared to the scale of the 't Hooft coupling of the Yang- 
Mills theory, the background becomes weakly curved everywhere. This was, in fact, the 
first model of a confining theory constructed following the AdS/CFT correspondence. 
In the limit of large N^ with fixed Nf, the back- reaction of the D8-branes is small, and 
they can be treated as probes in this background. Sakai and Sugimoto noted that in this 
setup there is a very nice interpretation of chiral symmetry breaking^. At large radial 
positions (corresponding to the high-energy limit of the field theory) the D8-branes and 
anti-D8-branes are separated and one sees the full U{Nf)i x U{Nf)R chiral symmetry. 
Apriori, as we go into the interior of space, the D8-branes and the anti-D8-branes could 
remain separated, or they could join together into Nf continuous 8-branes, so that only 
a single (diagonal) U{Nf) symmetry is visible, corresponding to spontaneous chiral sym- 
metry breaking. In the background of |2] , the radial position has a minimal value (where 
the X4-circle shrinks to zero size) and the 8-branes have nowhere to end, and therefore the 
D8-branes have to smoothly join the anti-D8-branes and exhibit xSB. In ^7\ an explicit 
solution for the position of the 8-branes was constructed (in the probe approximation) 
which realizes this. 

Recently, it was realized in J2I3 that the scale of chiral symmetry breaking in this 
model is actually independent of the scale of confinement (the first depends mostly on 
L while the second depends mostly on R), and that chiral symmetry breaking persists 
even in the R ^ oo limit in which confinement disappears^. We will see this separation 
of scales also in our finite temperature analysis. 

In this paper we will study the model of Sakai and Sugimoto at finite temperature. In 
asymptotically free gauge theories with confinement and chiral symmetry breaking, both 
of these properties disappear at high temperatures as the effective coupling becomes 
small. In QCD there is no sharp phase transition distinguishing the low-temperature 
and high-temperature phases, since there are no sharp order parameters for confinement 
(due to the presence of dynamical quarks) and for the (approximate) chiral symmetry 
breaking. However, in a theory with exactly massless quarks one expects to have a sharp 
phase transition associated with the chiral symmetry restoration (characterized by the 

^A similar structure of probe branes in a confining background was previously introduced in |S] in 
the context of the model of ^ . 

^The independence of the two phenomena in the gravity approximation was first noted for an Abelian 
chiral symmetry in |1HJ . 



chiral condensate as an order parameter), and in the large Nc limit the deconfinement 
transition also becomes a sharp phase transition (since the quarks can be ignored in 
this limit, and the pure gauge theory at finite temperature has a Z^^ global symmetry 
which is unbroken in the confined phase and broken in the deconfined phase). It is 
interesting to ask whether these two phase transitions occur at the same temperature 
or not. In the large A*";, limit of QCD this is a difficult question to answer analytically. 
However, in the Sakai-Sugimoto model (which is continuously related to large Nc QCD 
by varying parameters) we will see that the question is easy to answer^. We will find that 
simple geometrical arguments imply that chiral symmetry is always broken in the confined 
phase, so that chiral symmetry restoration can happen either at the same temperature as 
deconfinement or above it. A more detailed analysis will be necessary to decide between 
these two possibilities, and we will find that in the weakly-curved limit the answer to this 
question depends on the ratio L/R (which can go between zero and vr; in the limit of small 
L/R this ratio is inversely related to the ratio between the low-spin meson masses and the 
low-spin glueball masses). We will show that for L/R > 0.97 the two phase transitions 
occur together, while for L/R < 0.97 the chiral symmetry restoration transition happens 
at a higher temperature than the deconfinement transition. In the supergravity limit all 
of these phase transitions are of first order. 

We begin in section 2 with a description of the system at zero temperature. In par- 
ticular we present a short review of the Sakai-Sugimoto model and its meson spectrum, 
generalizing the construction of |^ to arbitrary values of L/R. In section 3 we discuss 
the finite temperature behavior of this theory. We start by discussing the deconfine- 
ment transition, which is associated with the bulk physics, and its implications on chiral 
symmetry breaking. We continue by studying the realization of chiral symmetry in the 
deconfined phase and computing the chiral symmetry restoration temperature. We also 
discuss the spectrum of quark- ant i- quark excitations of the theory at finite temperature. 
We end in section 4 with our conclusions. An appendix contains the computation of the 
difference of free energies which is relevant for the deconfinement transition. 

2 The zero temperature string theory model 

The string theory dual of the gauge theory we described above at zero temperature is 
the Sakai-Sugimoto model J7]. We start by reviewing the basic structure of this model. 
We then describe the spectrum of mesons of low and high spins. 



^Note that this theory is not asymptotically free, so apriori it is not clear whether at high temperatures 
it should deconfine and restore chiral symmetry or not. 



2.1 Review of the Sakai-Sugimoto model 

The model of Sakai and Sugimoto [1^ describes the near-horizon hmit of a configura- 
tion of Nc D4-branes, wrapping a circle in the X4 direction with anti-periodic boundary 
conditions for the fermions, and intersected by Nf D8-branes (located at X4 = 0) and 
Nf anti-D8-branes (located at x^ = L). This is dual to a 4 + 1 dimensional maximally 
supersymmetric SU{Nc) Yang-Mills theory (with coupling constant g^ and with a specific 
UV completion that will not be important for us), compactified on a circle of radius R 
with anti-periodic boundary conditions for the fermions, with Nf left-handed quarks lo- 
cated at X4 = and Nf right-handed quarks located at X4 = L (obviously we can assume 
L < TiR). 

In the limit Nf <^ A''^, this background can be described by Nf probe D8-branes 
inserted into the near-horizon limit of a set of Nc D4-branes compactified on a circle with 
anti-periodic boundary conditions for the fermions j2] . This background is simply related 
to the (near-horizon limit of the) background of near-extremal D4-branes by exchanging 
the roles (and signatures) of the time direction and of one of the spatial directions. Let us 
now briefly review this model, emphasizing the manifestations of confinement and chiral 
symmetry breaking. The background of type IIA string theory is characterized by the 
metric, the RR four-form and a dilaton given by 
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where t is the time direction and x* {i = 1,2,3) are the uncompactified world-volume 
coordinates of the D4 branes, X4 is a compactified direction of the D4-brane world- 
volume which is transverse to the probe D8 branes, the volume of the unit four-sphere 
^4 is denoted by V4 and the corresponding volume form by 64, Ig is the string length 
and finally gg is a parameter related to the string coupling. The submanifold of the 
background spanned by X4 and u has the topology of a cigar (as on the left-hand side of 
figure 121 below) where the minimum value of u at the tip of the cigar is u\. The tip of 
the cigar is non-singular if and only if the periodicity of X4 is 




5x4 = ^ P^ =27rR (2) 



and we identify this with the periodicity of the circle that the 4 -|- 1-dimensional gauge 
theory lives on. 

The parameters of this gauge theory, the five- dimensional gauge coupling g^, the low- 
energy four- dimensional gauge coupling 5(4, the glueball mass scale Mgi,, and the string 



tension Tgt are determined from the background (^ in the following form ^U] '■ 

T -J-/ 1 1 f uaY/^ _ 2 glN,_ As 

^* 27rZ2^^"^""'"="^ 27fP^ \RnJ 27nR'^ 277r2i?3' ^"^^ 

where A5 = S'fA^c, ^g6 is the typical scale of the glueball masses computed from the 
spectrum of excitations around (^, and Tgt is the confining string tension in this model 
(given by the tension of a fundamental string stretched at m = ua where its energy 
is minimized). The gravity approximation is valid whenever A5 ^ i?, otherwise the 
curvature at m ~ ma becomes large. Note that as usual in gravity approximations of 
confining gauge theories, the string tension is much larger than the glueball mass scale in 
this limit. At very large values of m the dilaton becomes large, but this happens at values 
which are of order N^'^ (in the large Nc limit with fixed A5), so this will play no role 
in the large Nc limit that we will be interested in. The Wilson line of this gauge theory 
(before putting in the D8-branes) admits an area law behavior 21j, as can be easily seen 
using the conditions for confinement of P2] . 

Naively, at energies lower than the Kaluza-Klein scale 1/ R the dual gauge theory is 
effectively four dimensional; however, it turns out that the theory confines and develops 
a mass gap of order Mgi, = 1/R, so (in the regime where the gravity approximation is 
valid) there is no real separation between the confined four-dimensional fields and the 
higher Kaluza-Klein modes on the circle. As discussed in J2], in the opposite limit of 
A5 ^ R, the theory approaches the 3 + 1 dimensional pure Yang-Mills theory at energies 
small compared to 1/R, since in this limit the scale of the mass gap is exponentially small 
compared to 1/R. It is believed that there is no phase transition when varying A5/-R 
between the gravity regime and the pure Yang-Mills regime, but it is not clear how to 
check this. 

Next, we introduce the probe 8-branes which span the coordinates t, a;*, ^4, and follow 
some curve u^x^) in the (0:4, u)-plane. Near the boundary at u —^ 00 we want to have 
Nf D8-branes localized at X4 = and Nf anti-D8-branes (or D8-branes with an opposite 
orientation) localized at X4 = L. Naively one might think that the D8-branes and anti- 
D8-branes would go into the interior of the space and stay disconnected; however, these 
8-branes do not have anywhere to end in the background (Q), so the form of u{x4) must 
be such that the D8-branes smoothly connect to the anti-D8-branes (namely, u must 
go to infinity at 0:4 = and at X4 = L, and du/dx^ must vanish at some minimal u 
coordinate u = uq). Such a configuration spontaneously breaks the chiral symmetry 
from the symmetry group which is visible at large u, U{Nf)i x U{Nf)fi, to the diagonal 
U{Nf) symmetry. Thus, in this configuration the topology forces a breaking of the chiral 



symmetry; this is not too surprising since chiral symmetry breaking at large Nc follows 
from rather simple considerations [231 • 

In order to find the 8-brane configuration, we need the induced metric on the D8- 
branes, which is 
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where u' = du/dx^. It is easy to check that the CS term in the D8-brane action does not 
affect the solution of the equations of motion. More precisely, the equation of motion 
of the gauge field has a classical solution of a vanishing gauge field, since the CS term 
includes terms of the form C^ A F A F and C3 A F A F A F. So, we are left only with 
the DBI action. Substituting the determinant of the induced metric and the dilaton into 
the DBI action, we obtain (ignoring the factor of A'^^ which multiplies all the D8-brane 
actions that we will write) : 



Sdbi = Ts dtd^xdxid'^Qe '^ y — det{g) = — / dx^u"^^ 
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where g is the induced metric (JH) and Tg includes the outcome of the integration over all 
the coordinates apart from dx^. 

The simplest way to solve the equation of motion is by using the Hamiltonian of the 
action (jS)), which is conserved (independent of X4) : 
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where on the right-hand side of the equation we assumed that there is a point Uq where 
the profile u{x^ of the brane has a minimum, u'{u = uq) = 0^. We need to find the 
solution in which as u goes to infinity, X4 goes to the values X4 = 0, L; this implies 



dxA 
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with u' given (as a function of u) by (jH)) (note that m is a double-valued function of X4 in 
these configurations, leading to the factor of two in ((Tj)). The form of this profile of the 
D8 brane is drawn in figure Q^a). Plugging in the value of u' from (jHI) we find 

'■°° du „q/9 r°° . 1 
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®This type of analysis was done previously for Wilson line configurations. See, for instance, I21I. 
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Figure 1: The dominant configurations of tlie D8 and anti-D8 probe branes in tlie Sakai- 
Sugimoto model at zero temperature, wliicli break tlie cliiral symmetry. The same con- 
figurations will turn out to be relevant also at low temperatures. On the left a generic 
configuration with an asymptotic separation of L, that stretches down to a minimum at 
u = uq, is drawn. The figure on the right describes the limiting antipodal case L = ttR, 
where the branes connect at Uq = u\. 
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where |/a = ua/uq. Small values of L correspond to large values of Uq. In this limit we 



have y\ <^ 1 leading to L oc JR^^/uq. For general values of L the dependence of Uq on 
L is more complicated. 

There is a simple special case of the above solutions, which occurs when L = nR, 
namely the D8-branes and anti-D8-branes lie at antipodal points of the circle. In this 
case the solution for the branes is simply X4{u) = and X4{u) = L = ttR, with the two 
branches meeting smoothly at the minimal value u = uq = u^ to join the D8-branes 
and the anti-D8-branes together. This type of antipodal solution is drawn in figure 
CJb). It was shown in \T2\ that this classical configuration is stable, by analyzing small 
fluctuations around this configuration and checking that the energy density associated 
with them is non-negative. 

For general L, one might apriori think that there could be two ways for the D8- 
branes to connect to the anti-D8-branes - one in which they connect by curving towards 
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positive X4 (in the direction where the D8-anti-D8 distance is L) and the other by curving 
towards negative X4 (in the direction where the distance is 2ttR — L). However, smooth 
configurations only occur when the integral (jTj) is less than or equal to ttR, so there is 
only one smooth configuration for every value of L. 

Note that the Sakai-Sugimoto model has 3 dimensionful parameters : A5, L and R, 
and gravity is reliable whenever A5 ^ R. The physics depends on the two dimensionless 
ratios of these two parameters; we will see that in the gravity limit the mass of the (low- 
spin) mesons is related to 1/L while the mass of the (low-spin) glueballs is related to 
1/R. As discussed above, in the limit A5 <^ i? this theory approaches (large Nc) QCD at 
low energies. This remains true also after adding the fiavors, at least when L is of order 
R (note that the parameter L is not visible in the four-dimensional low-energy effective 
action). However, this QCD limit is opposite from the supergravity limit, and there are 
no known reliable methods to analyze it. 

2.2 The spectrum of mesons 

The mesons of the Sakai-Sugimoto model are described by strings ending on the D8- 
branes. It is convenient to separate the discussion to low-spin mesons and high-spin 
mesons. The former are described via modes of the massless fields living on the D8- 
branes, whereas the latter are associated with string configurations that "fall" from the 
D8-branes down to the wall at u = ua, stretch along the wall and then go back up again 



Let us start with the low-spin mesons. These mesons correspond on the string theory 
side to the fiuctuations of the massless fields on the probe branes. The fluctuations of 
the gauge fields on the branes give pseudo-vector and scalar mesons as well as massless 
"pions" , and the fiuctuations of the scalar field describing the embedding of the branes 
give massive scalar mesons. In P7| these two sectors of the mesonic spectrum were 
computed for the special case of mq = ^^a, namely L = nR. It is not difficult to generalize 
this analysis to general values of L which we are interested in here, using the analysis of 
the fiuctuations performed in ^3j which we summarize below (the discussion in ^3] is 
for a non-critical AdS^ model, but one can easily apply the same analysis for our model). 

As an example we describe in detail the modes coming from the components of the 
gauge field living on the D8-branes which are not in the S"^ directions. If we write the 
action using the coordinate u (instead of using X4 as in (0)), the equations of motion for 



these fields that follow from the DBI action are given by^ 



u'/^^iuy/^d^F'^'' -duiu'/^liuY^^F"^) = 0, (9) 

d^F^"" = 0, (10) 



where /x, t/ = 0, 1, 2, 3 and where we have defined 



u 



We now separate variables and expand the five dimensional gauge fields in modes in the 
radial direction as follows : 

n n 

Using this decomposition, equation (|Tn|l reads: 

Y: {my"\xn<PM - {d,B^^^{xn)duM^^)) = O, (is) 

n 

where m„ is the four-dimensional mass of the field vr*^"'', namely r]^'^d^d^7r^^^ = m'^7i^"'\ 
We can now choose a Lorentz gauge c^^A^ = — >■ d^B'fs = 0, and in this gauge it is easy 
to see that all the n^"''' fields are set to zero by (fT^ except for the Goldstone bosons tt^^^ 
with rho = 0, which survive as a massless field for any value of Uq. The other equation 
of motion Q fixes the form of (pQ to be 

00 oc ^7^ . (14) 



u 



5/2 ■ 



Finally, the eigenvalue equation for the modes ipn{u) reads 

where rrin are the masses of the four- dimensional vector fields 5*^"-' , such that their equa- 
tion of motion is 'q'^^d^.dpBJj!^^ = rn^B'^J^\ Imposing appropriate orthonormality conditions 
[Ej , one finds that the effective four-dimensional action for these fluctuations of the gauge 
field living on the probe D8-branes can be written as: 



S = — d xtr 






(16) 



^We perform the analysis for a single D8-brane, but the generalization to Nf D8-branes is simply 
performed by turning all the fields into Nf x Nf matrices, turning the derivatives into [/(A^^)-covariant 
derivatives, and taking a trace in the action. 
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The fields vr'-^'' are the Goldstone bosons associated with the spontaneous breaking of the 
U{Nf)L X U{Nf)ji global symmetry to the diagonal U{Nf) (which is the global symmetry 
of the remaining action ()16p). As in ^Tj, one can also go to the Au = gauge, where 
there are no vr'^") fields and the Goldstone bosons are encoded in the boundary value (as 
u — > oo) of y4^. Note that the pions (the n^^^ fields) are massless for any value of uq (or 
for any value of L), as expected since the quarks are massless (in the sense of having a 
vanishing current algebra mass) for any value of L. 

The explicit mass spectrum of the vector mesons is found by looking for normalizable 
eigenfunctions of (fTHjl . In the special case of uq = u\ this spectrum was derived in pTj . 
For generic probe profiles, for which uq ^ u\, the spectrum was computed in jHj for a 
similar non-critical string model. In general this spectrum is quite complicated. In the 
limit of Mq ^ Ma or L <^ i? the spectrum simplifies. In this limit we can approximate 
f{u) ^ 1, and there is a scaling symmetry of (fTH|) which implies that all m^ are pro- 
portional to uo/-Ri)4. Equation (jHl) then implies that this is proportional to 1/L^, so in 
this limit all the low-spin meson masses are simply proportional to 1/L. This illustrates 
the separation of scales between the low-spin mesons and the low-spin glueballs (whose 
masses are proportional to 1/R in the supergravity approximation) in the Sakai-Sugimoto 
model (as emphasized in pUj). 

The DBI action of the D8-branes only includes mesons of spin one or less. To describe 
higher-spin mesons we need to look at more general string configurations that start and 
end on the probe branes. For large spin these strings are long and can be described semi- 
classically. The relevant string configurations can be schematically decomposed into three 
parts : a segment from the probe brane (at u = uq) to the wall {u = u\), then a segment 
that stretches along the wall in the spatial directions, and then another "vertical" part 
stretching from the wall back to the probe brane, as depicted in figure |21 

By computing the energy of the string it was shown in J2] that this string is equivalent 
to a string in fiat space-time with two massive endpoints. The mass of the end-points, 
which can be identified with the constituent quark mass m^ , is given by [T^ 



j_ /■"() 1 /•«() 



which is simply the mass of a string stretched from the D8-branes to the minimal value of 
u. Note that for small L this constituent mass scales as uq/o' oc As/L^. Let us emphasize 
again that the bare mass and current algebra masses are always zero in our configurations, 
and there is no simple description of the deformation associated with turning on a bare 
mass for the quarks (which corresponds to turning on a non-normalizable mode of the 
string connecting the D8-branes to the anti-D8-branes). Using the map to the string 
with massive endpoints it is straightforward to calculate the classical mass and angular 
momentum of these stringy mesons ^21 121| • Somewhat surprisingly, it was found in 
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Figure 2: The high-spin "stringy" meson is a string starting at the lowest point of the 
probe branes u = uq, going down to the wall, stretching horizontally in space along the 
wall, and then going back up vertically to the probe branes at u = Uq. 

(in a similar model to the one we describe) that for mq of order u\ there is a range of 
values of Uq/ua for which the masses of the light vector and pseudo scalar mesons scale 
roughly linearly with m^ , reinforcing its identification as a constituent quark mass. The 
decay processes of the stringy mesons were discussed in [T3]. 

3 Thermodynamics of the Sakai-Sugimoto model 

3.1 Review of the bulk thermodynamics 

In the large N^ limit (with finite Nf), the effect of the D8-branes on the background is 
subleading; the difference in free energies between different bulk backgrounds is of order 
A^^, while the contribution of the D8-branes is of order Nc * Nf. Thus, we can analyze 
the bulk thermodynamics independently of the presence of the D8-branes, and then add 
them in as probes to the dominant bulk background at each temperature. 

The thermo dynamical behavior of the bulk is analogous to the thermodynamical be- 
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Figure 3: The topology of the solutions which dominate in the low temperature (confined) 
and high temperature (deconfined) phases of the Sakai-Sugimoto model, as refiected in 
the {t, u) and (0:4, u) submanifolds. 

havior of other theories compactified on circles with anti-periodic boundary conditions for 
fermions in the gravity approximation (see, e.g., |21I2H1)- In the gravity approximation, 
in order to study the theory at a finite temperature T, we need to look at Euclidean grav- 
itational backgrounds which are asymptotically equal to (0), but with the time direction 
Euclidean and identified with a periodicity /3 = 1/T, and with anti-periodic boundary 
conditions for the fermions along this direction (as well as along the X4 direction). One 
such background is obviously given by the Euclidean continuation of (^ itself, with t 
periodically identified; it is easy to verify that this is a smooth background which is 
weakly curved (at least until the temperature reaches a number of order the square root 
of the confining string tension T^t ~ X^/R^, at which stage the size of the time circle at 
the minimal radial position u = u\ becomes of order the fundamental string scale and a 
tachyon develops). In this background the periodicity of t is arbitrary and equal to /3, 
and the periodicity of X4 is 2'kR with R related to ma by (j21). 

In the background described above the X4 circle shrinks when u reaches its minimal 
value, but the t circle never shrinks to zero size. Another background with the same 
asymptotics is given by exchanging the behavior of the t and x^ circles, or equivalently 
by moving the factor of f{u) in (P) from the dx^ term to the dt^ term. Now, the t 
circle shrinks to zero size at the minimal value of u (which is now related to (3 rather 
than to 277 i?, and given by ut = 167r^i?|,4/9/?^), while the x^ circle never shrinks. This 
background also exists as a classical solution for any value of the periodicities of X4 and 
t (and develops a tachyon for very small values of R). 

These two backgrounds are the only known smooth Euclidean solutions which asymp- 
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totically approach ^ ^°. Thus, in order to decide which background dominates at a given 
temperature T we need to compute the free energies of these backgrounds, given (in the 
gravitational approximation) by the classical supergravity action times the temperature, 
and see which one has a lower free energy. As usual, the classical action actually di- 
verges and needs to be regulated. One way to regulate it is by looking at the difference 
between the actions of the two solutions, which turns out to be finite (this method is 
not always reliable, but it gives the correct answer in this case). Another way is to add 
counter-terms at a large u cutoff that make the action finite as the cutoff goes to infin- 
ity. In the first method it is clear that the free energies of the two solutions agree when 
the (asymptotic) circumferences of the two circles are equal, (3 = 27iR, since the two 
solutions are identical then except for a relabeling of the coordinates. Thus, at this tem- 
perature Trf = 1/2'kR there is a first order phase transition between the two backgrounds 
(the transition is of first order since the solutions do not smoothly connect there, but 
continue to exist as separate solutions both below and above the transition). It is easy 
to see that the background in which the 0:4 circle shrinks to zero size dominates at low 
temperatures T < 1/2'n'R, while the background in which the t circle shrinks dominates 
at high temperatures, T > 1/2ttR; intuitively, the circle which is smaller prefers to shrink 
to zero size. The relevant computation is described in the appendix, and leads (jHEl) to a 
difference of free energy densities proportional to N^ R{gsNc)[{2TTT)^ — 1/R^]- 

The physical interpretation of this phase transition (which is a close relative of the 
Hawking-Page transition in asymptotically AdS backgrounds 0) is straightforward. If 
we compute the quark-anti-quark potential in the two backgrounds (before adding the 
D8-branes) using the methods of [SEl I2Z1 122] , we find that in the low-temperature back- 
ground in which y/gugxx is finite at the minimal value of u the potential is linear in 
the distance, corresponding to confinement, while in the high-temperature background 
in which ^JgttQxx goes to zero at the minimal value of u the potential decays with the 
distance, corresponding to a deconfined phase. Similarly, a computation of the value 
of the free energy in the two phases (which requires adding appropriate counter-terms) 
yields a result of order N^ in the low-temperature phase, and a result of order A^^ in the 
high-temperature phase. 

As mentioned above, as we increase R compared to the scale set by the five dimen- 
sional gauge coupling A5, the supergravity background becomes highly curved, and the 
theory at low energies approaches the four dimensional pure Yang-Mills theory. It is be- 
lieved (though it is difficult to prove) that the deconfinement transition described above 

^°Note that asymptotically these backgrounds become strongly coupled so we should really lift them 
to M theory, and discuss the asymptotics in the M theory regime, but in the large Nc limit the strong 
coupling region can be arbitrarily far away from the regions of finite u which we are interested in so this 
does not play any role. 
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Figure 4: A minimal conjecture for the schematic form of the phase diagram of the 4 + 1 
dimensional maximally supersymmetric SU{Nc) gauge theory compactified on a circle 
with anti-periodic boundary conditions for the fermions, as a function of the dimensionless 
coupling constant X5/R and the dimensionless temperature TR, based on the known 
results in the various limits and on the T ^->- 1/2ttR symmetry. All solid lines in the 
diagram denote phase transitions; the dashed line may or may not be a phase transition 
line. 



is connected to the deconfinement transition of the large Nc pure Yang-Mills theory in 
this limit. However, this connection is somewhat subtle [2Hl! since even for small \^/R it 
is clear by the symmetry arguments described above that there is always either a phase 
transition line at T = 1/2ttR, or else the dominant phase at this temperature must be 
invariant under the exchange T ^^ 1/2ttR (which is not the case for the confined and 
deconfined phases discussed above). On the other hand, the deconfinement phase transi- 
tion in the pure Yang-Mills theory occurs at a temperature of order Aqcd ^ 1/R. Due 
to the T ^^ 1/2ttR symmetry there must then be at least one more phase transition line 
at a temperature much larger than 1/R. Thus, we see that new phases must appear as 
we approach the four dimensional Yang-Mills limit. The phase structure of this theory 
in that limit includes several phase transition lines, which, presumably, all join into the 
line T = l/27ri? in the gravity limit. A minimal conjecture for the phase structure is 
depicted in figure |31 note that with this conjecture the deconfinement transition line 
in the supergravity limit is connected to the deconfinement transition line in the four 
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dimensional Yang-Mills limit, but this connection does not seem to be smootli^^. 

3.2 The low temperature phase 

As described above, the background corresponding to the low-temperature phase of the 
theory is identical to that of the Sakai-Sugimoto model at zero temperature ((H), apart 
from the fact that the time direction is Euclidean and compactified with a circumference 
P = 1/T, where T is the temperature of the dual gauge theory. Hence, just as for the zero 
temperature case, the dual gauge theory is in the confining phase and the string tension 
is given by Q. Furthermore, the solution for the profile of the D8 brane is still given 
by (jHI)- Note that this agrees with general arguments that the chiral condensate (in the 
large Nc limit) should be independent of the temperature in the confined phase JHIIIIS]- 
Defining y = u/uq, y^ = ua/uq, and f{y) = 1 - {yK/yf (such that /(I) = 1 - |/i), we 
find from (jH)) that 



u \ 



3/2 



Substituting this outcome of the equation of motion, the action is now 



J{y) 



f,C'3/2 7/2 5/2 



2f.R''luV^ '•I 



1-04 "0 



v^yWi? 



dz ■ (19) 



3^, Jo ^13/6 /l _ (1 _ ^3)^8/3 _ ^3^ 



The form of the profile ()18|) of the D8-branes is drawn in figure ^a). Again, the 
topology of the background forces a profile that smoothly connects the D8-branes to the 
anti-D8 branes, signaling chiral symmetry breaking in the dual gauge theory. 

The relation between L and uq is precisely the same as for zero temperature (jH)). For 
small values of L we find that the dependence of the action (|19|). which is proportional 
to the vacuum energy of the system, on the separation distance is given by 

>3db/ oc —^. (20) 

For general values of L the dependence is more complicated. 



^^This conjectured phase diagram is the same as the phase diagram of the same theory reduced on the 
three spatial directions which was discussed in [221, suggesting that perhaps the phase structure does 
not change when the three infinite spatial dimensions are compactified on a torus. We thank M. Van 
Raamsdonk for discussions of this point. 
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3.3 The high temperature phase 

The gravity background associated with the high temperature phase of the dual gauge 
theory takes the form 



ds' 



^(4) 



U 



2'KNr 



3/2 



f{u)d£' + dijdx'dx^ + dx 



+ 



R 



Di 



U 



3/2 



u'^dQl + 



du" 



fiu) 



V, 



-€4, 



u 



R 



D4 



3/4 



Rl, = ng^NJl f{u) = 1 



^V>(21) 



u J 



in the same notations as in (^). The time direction t now shrinks to zero size at the 
minimal value of m, m = mt, and in order to avoid a singularity there the time direction 
must have a period of 



* = T 



4i IR'; ^ "^ 



•-04 



Ur 



/5 



(22) 



which we identify with the inverse temperature of the system /? = 1/T. On the other 
hand, in this solution the x^ direction can have an arbitrary periodicity I'kR. 

We expect that again the D8 probe branes span the coordinates t,x\r24, and are 
described by some curve u{x4). The induced metric on the D8 branes now reads 






+ 



u \3/2 
Roi/ 



R \ 3/2 



U 



f{u)dt^ + Sijdx'dx^ 
u^dVl 



+ 



u 



R 



DA 



3/2 



1 + 



Rda\^ m'^ 



u 



') fiu) 



(XJb A 



(23) 



where u' = du/dx^. Just as for the low temperature phase, here too the CS term does 
not affect the solution of the equations of motion. Upon substituting the determinant of 
the induced metric and the dilaton, the DBI action is given by 



Sdbi = Ts dtd^xdxid^Vte' 



det(^) 



9s 



dx^u^Jfiu] 



\ 



1 + 



JT/)4 \ U 



U 



f{n 



, (24) 



where Tg includes the outcome of the integration over all the coordinates apart from dx^. 
Again, conservation of the Hamiltonian of (j24j) implies that 



f{u)u' 



^^ \ u ) f{u) 



constant = UqJ f{uo) 



(25) 



where on the right-hand side of the equation we assumed that there is a point Uq where 
the profile of the brane has a minimum, u'{u = mq) = 0. Defining again y = u/uo, we 
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Figure 5: The possible configurations of tlie D8 and anti-D8 probe branes in tlie liigli 
temperature (deconfined) pfiase. A generic connected configuration with an asymptotic 
separation of L, that stretches down to a minimum at m = mq, is drawn in (a). This 
corresponds to the chiral symmetry broken phase. Figure (b) depicts a chiral-invariant 
solution, with separate D8 and anti-D8-branes. 



find that 



u 



u 



R 



DA 



3/2 



f{y)^ 






2/8-1. 



Substituting the outcome of the equation of motion, the action is now 



S 



DBI 



9s 



dy 



,5/2 



1 



/(I) 



(26) 



(27) 



Since we assumed u'{u = uq) = 0, this action corresponds to a configuration where the 
D8-branes merge smoothly into the anti-D8-branes, as in figure El^a). 

As discussed above, this configuration exhibits in the string theory dual the phe- 
nomenon of xSB in the gauge theory. However, there is yet another classical config- 
uration that solves the equations of motion. This is a configuration of a constant x^, 
X4{u) = and X4{u) = L (or u' -^ oo), see figure Efb). This configuration is possible 
in the high-temperature phase since in this phase the t circle shrinks to zero aX u = ut, 
so the D8-brane can just smoothly end there (and wrap around the "cigar" in the t — u 
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Figure 6: AS" as a function of yx, in units of 2T^R/)^Uq / Qs- 

plane); in the Minkowski space continuation of this solution, the D8-brane would go into 
the horizon of the black D4-branes. For such a configuration that naturally describes the 
situation of an unbroken chiral symmetry, the action reads (summing over the D8 and 
anti-D8 contributions) 



r-O 
'-'DBI 



2T8 
9s 



duu 



ux 



R 



DA 



U 



3/2 



9s 



'dyy"'+ dyy'/' 



Vt 



(2^ 



where yr = ut/uq. 

We are now facing the question of which of the two possible configurations is preferred. 
This is determined by the free energy, or the action, of the configurations. The relevant 
quantity is the following difference between the two actions 



AS 



9s{S-S^) 
2fsR'jiyj' 



dyy'/' 



1 



i-IS^- 



dyy 



5/2 



Vt 




l-y^z 



^1 _ ^3^ _(1_ ^3)^8/3 



^^{i-vT) 



(29) 



We are unable to evaluate AS analytically so we do it numerically. The result is 
drawn in figure IHl which shows AS as a function of yx- It is easy to see from this figure 
that for yT > yx ^ 0.73572, AS' is positive and hence the preferred configuration is that 
of disconnected D8 and anti-D8-branes, whereas in the region where 2/t < l/r ~ 0.73572 
the system will be dominated by a smooth configuration that connects the D8 and anti 
D8 branes, namely the chiral symmetry broken phase. 
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Next, we would like to express the critical point in terms of physical quantities. Recall 
that yr = ut/uq. ut is related to the temperature in the gauge theory by (j221)- Uq is the 
minimal point of the connected probe brane configuration, which can be related to the 
asymptotic distance L between the positions of the D8 and anti-D8-branes (see figure Ej) 
in the following manner^^ 




(30) 

Note that for small L, namely for yx « 1, we find Uq oc R^^/L"^ and again S oc 1/L^. 
We would like to relate the asymptotic distance L at the transition to the temperature of 
the transition. Evaluating the integral in ()H()|1 at the transition temperature y^ ~ 0.73572 
gives Lc — 0.751(i?|,4/uo)^''^, and using equation (f^ we find 

T^SB ^ 0.154/L. (31) 

Recall that the deconfinement phase transition is at T^ = 1/2ttR, and that this 
analysis applies only in the deconfined phase. For L > 0.97R the temperature in the 
deconfined phase is always higher than (jHTj) . so the deconfinement and chiral symmetry 
restoration transitions happen together. On the other hand, for L < 0.97R, the two phase 
transitions are separate, and there is a finite range between Td and T^sb for which the 
theory is deconfined but with a broken chiral symmetry. The full phase diagram of the 
theory (in the supergravity approximation) is drawn in figure |7| Note that even though 
for the gravity approximation to be good everywhere in space we required A5 ^ R, for 
gravity to be good near the 8-branes it is actually enough to require A5 ^ L, and at high 
temperatures gravity gives a good approximation everywhere whenever A5 ^ 1/T. Thus, 
we can trust our analysis of the phase transition whenever A5 ^ L. In particular, in the 
R ^ 00 limit discussed in [20j, the theory is always deconfined, and passes through a 
chiral symmetry restoration transition at T = T^sb given by (jHT|) . ^^ 

3.4 The spectrum in the various phases 

In this subsection we discuss the spectrum of quark configurations in the different phases 
described above; such a discussion is more appropriate in Minkowski space, in the 

^^In fact, as we discussed above, this integral could equal either L or 2ttR — L, but it is easy to see 
that the second option is always sub-dominant wherever it exists. 

""^^Note added : in parallel with this paper, the paper [32] appeared which discusses the same system 
in the R —^ 00 limit, with identical results to ours. 
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Figure 7: The phase diagram of the Sakai-Sugimoto model at finite temperature, in the 
gravity approximation A5 3> R. In this approximation the phase structure depends only 
on the two dimensionless parameters TR and L/ R. For small L/ R the deconfinement 
and chiral symmetry restoration transitions are separate, while for L/ R > 0.97 there is 
a single transition. 

finite-energy-density configurations (in the microcanonical ensemble) obtained by Wick- 
rotating the Euclidean configurations discussed above. In general there are two types of 
possible configurations. We can have strings with both ends on D8-branes, corresponding 
to mesons (for light strings these are described by the DBI action, and for large spin they 
can be described by semi-classical strings, see figure 12)). And, if the background contains 
a horizon (as in the deconfined phase) we can have a string starting on an 8-brane and 
going into the horizon, corresponding to a deconfined quark (or anti-quark). 

As discussed above, the behavior in the low-temperature phase is very similar to the 
behavior at zero temperature, as expected in a theory with a small number of weakly 
interacting massless particles (the "pions"). There is a discrete spectrum of mesons, and 
there are no free quark states. Computing the spectrum of mesons (including the pseudo 
vector and scalar mesons described in §2.2) and glueballs in this phase gives exactly the 
same results as for zero temperature. The high spin mesons are described as classical 
string configurations stretched from the D8 probe brane down to the "wall" at u = u\ 
and back to the probe branes P^ (see figure |2I), and again there is no difference between 
the zero temperature and low temperature results. 

As we raise the temperature, at T = 1/2ttR we go over to the deconfined bulk 
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phase; as we saw, depending on L/ R, we can either get a phase with broken chiral 
symmetry or with unbroken chiral symmetry. In the case of small L/ R with broken 
chiral symmetry the D8-branes are still connected. In this phase the spectrum of bulk 
modes has no mass gap, due to the presence of the horizon in (J2I])- On the other hand, 
a computation of the meson spectrum using the D8-brane action still gives a discrete 
spectrum as in [TT]. However, unlike the low-temperature case, here the spectrum of 
quark- ant i- quark configurations is continuous at high enough energies. This is because 
this spectrum includes strings starting on the D8-brane and going down to the horizon 
u = ut (in Minkowski space such configurations are just deconfined quarks, while in 
Euclidean space these strings must go back to the D8-brane, but they can stretch any 
distance in the three dimensional space before doing so with no cost in energy). These 
excitations give rise to a continuous spectrum of quark-anti-quark configurations, which 
starts above a mass gap which equals twice the mass of a string stretching from uq to ut-. 



Mgap = 2J^° ^1 - {uT/ufdu/ll. This is obviously the analog of 2m^ defined in (fTTfl . 
and this also scales as As/L^ for small L/ R. Note that after we put in the D8-branes we 
can no longer use the spatial Wilson loop as a diagnostic for confinement, since this can 
be screened by the dynamical quarks, but the scaling of the free energy at large N^. can 
still be used to distinguish the confining phase at T < 1/2'kR from the deconfined phase 
at T > 1/2ttR. 

Even when L/ R is small and such an intermediate phase exists, at high enough 
temperature we always end up in a deconfined phase with a restored chiral symmetry, in 
which the 8-branes end at the horizon u = ut- In this phase there is no longer a mass gap 
for quark-anti-quark configurations, since strings stretched between the 8-branes and the 
horizon can be as light as we want, and the spectrum of quark-anti-quark configurations 
is continuous (just like the bulk spectrum). 

4 Summary 

In this paper we computed the phase structure of the Sakai-Sugimoto model, which is 
an example of a confining gauge theory with chiral symmetry breaking, related by a 
continuous change of parameters to large N^. QCD. Our results are summarized in figure 
[7| We found that the topology of the solutions (in the supergravity approximation) 
implies that the confined phase must always break the chiral symmetry. On the other 
hand, in the deconfined phase a computation is required to see whether chiral symmetry 
is restored at the deconfinement temperature or above it, and we found that in some 
range of parameters the first possibility is realized, while in another range the second 
possibility occurs. When L is much smaller than R we found that the chiral symmetry 
restoration scale is much higher than the deconfinement scale, as expected because of 
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the separation between the confinement and chiral symmetry breaking scales in this case 
(as emphasized in |2D1)- In this case the chiral symmetry restoration happens at a scale 
related to the mass scale of the mesons l/L, while the deconfinement transition happens 
at a scale related to the glueball mass scale 1/R. 

The fact that confinement implies chiral symmetry breaking is true also in a larger 
class of holographic duals of confining theories. When the confinement is realized geo- 
metrically (this is not always true jS]) there is a minimal value of the radial coordinate 
u = ua, and some spatial cycle shrinks there. A chiral symmetry should be realized by 
having two sets of branes (each carrying one U{Nf) factor), and chiral symmetry break- 
ing should be realized by the branes becoming connected. If the branes can end without 
reconnecting, confinement does not imply xSB; clearly this can only happen if the branes 
are wrapped on the cycle which vanishes at the minimal value u = u^. Conversely, if the 
branes are not wrapped around this cycle (as is the case in the Sakai-Sugimoto model 
which we discussed here), they do not have anywhere to end, so the chiral symmetry 
must be broken in the confined phase. It would be interesting to relate this discussion 
to the gauge theory discussions of chiral symmetry breaking in large N^ gauge theories 
[HOI IM] which show that chiral symmetry must be broken in the confined phase in the 
large Nc limit. 

In this paper we discussed the phase structure of the theory only in the gravity regime 
A5 ^ L. It was recently realized in [201 that there is another regime in which the theory 
is under control, which is the weak-coupling regime A5 <C L with R going to infinity. 
It would be interesting to study the phase structure also in this regime. There are also 
many other possible generalizations of the Sakai-Sugimoto model which can be analyzed 
by similar methods to the ones we used here, such as having several different stacks of 
D8-branes (and of anti-D8-branes) separated in the X4 direction. 

It is interesting to compare our results with the situation in large Nc QCD, which is 
related to our theory by taking the limit of small A5/-R. The technology to study large Nc 
QCD with chiral fermions on the lattice has only recently been developed, and preliminary 
results presented in jS^ suggest that in this theory chiral symmetry restoration occurs 
at the same temperature as the deconfinement transition, but that if one supercools 
the deconfined phase below the deconfinement temperature one encounters a critical 
temperature where the chiral symmetry is broken. This is similar to what we find in 
the case of L/R > 0.97 (where we also find a transition in the supercooled deconfined 
phase at the temperature Tc = 0.154/L), so it is plausible that the phase diagram in this 
regime connects continuously to that of large Nc QCD (up to the caveats discussed at 
the end of ^3.1). 
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A Appendix: The bulk free energies of the low and 
high temperature phases 

In this appendix we present the calculation of the difference between the free energies of 
the two possible solutions for the finite temperature description of the bulk, described in 
section 3.1. A similar calculation in the context of the AdS^ x S^ model was done in [2j. 
Rather than computing the action directly, we translate our results into the notations 
of [S3 and use the results for the action computed there. Thus, we parameterize the 
Euclidean metrics (0) and ()21|) that we are interested in as 

Ifds-" = dr^ + e^^^^Uxl + e^'^^^Uxl + e^'^^^d^l (32) 

where Xc is either x^ or the Euclidean time (the one whose circle shrinks to zero size 
at the minimal value of m), x\\ are the other four coordinates of i?^'^ (one of which is 
also compactified), r is the radial direction and dVL\ is the metric on the transverse 5*^. 
The background includes in addition the dilaton and a four-form RR field. For such 
backgrounds that depend only on a radial direction, the full type IIA supergravity action 
reduces to the following + 1-dimensional action 



S = V j dp (-4(A')' - (A')' - 4(z/')' + {^'f + 126-''^-'^ - g2e4A+A-4.-^j ^ ^33^ 



where the new radial coordinate p is defined by dr = —e^'^dp, primes denote derivatives 
with respect to p, and (f = 2(j) — AX — X — 4z/. V is the volume of all other directions 
(except p) in string units; this includes the volume V3 in the three spatial directions, 
the volume of the two circles which is 27ri?/5//^ and the volume of a unit 5*^ which is 
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27r^. The last term in the action corresponds to the contribution of the RR flux, with 
Qc = SttNc/V^ (see, for instance, jTHI) where Nc is the quantized RR flux. 
The equations of motion associated with the action ()33p are 



A" _ ig2g2(4A+A- 
A" _ lQ2g2(4A+A- 



2„2(4A+A- 



U 



':Qie 



2g2(4A+A+3i/-2</.) I j|:Q2g2(4A+A-</-) 



0, 
0, 
0, 
0. 



(34) 



It is straightforward to check that the following are the solutions of these equations 
, EHl corresponding to the background (Hj) : 



A 



-iln(sinh(-4M) + ^ln(^; 



A + ln{gs 



bp, 



X = X + Abp 
A 



u 



+ \n{Rn4/ls), (35) 



where b = —Q cu\/ A\/2R\)^g g. Here we wrote the solution of the low temperature phase, 
for the high temperature phase we have to replace ma with ut- 

Since the backgrounds (^ and ^Z\\ are expressed in terms of the radial coordinate 
u, it is useful to rewrite the action ()33|) in terms of integrals over m, which run from u^ 
to oo for the low temperature phase and from ut to oo for the high temperature phase. 
The action ()33|) then reads (with dots denoting derivatives with respect to u) 



S 



-V / du 

JuA 



-4A2 - X 



Au'^ + if"' 



du 
dp 



Vie 



8A+2A+6;/- 



Q\e 



2^8A+2A-20 
c 



dp 
du 



where the minus sign arises because dp/du is negative. Substituting the solution 
action takes the form 



(36) 
, the 



S 



-V / du 

JuA 



— Il 1 



du 
dp 



12 



Ql 



.gtn 9lR 



u'f{u) 



s -" 1-04 / 



|M37) 



At this point we need to compute the derivative -r- and its inverse. Using the relation 



g2A/g2A = [1 _ (!£A)3] ^nd ^ wc find that ^ is given by 

3 1 



dp 
du 



8^«[(^)'-l]' 



(3^ 
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Substituting this expression into the action we find that the divergence at large u is 
independent of ua, so it makes sense to subtract the expressions with u\ and with ut 
and to obtain a finite answer. The result for the difference between the action densities 
in the low temperature phase and in the high temperature phase is given by (defining 
b = bl^/u\ which is a constant independent of u/Cj 



AS S, 



low 



s, 



high 



V, 



V, 



An^Rp 5 



3 46 [gt 2g^sRD,. 



(ul - 4) 



(4 - o 



A. 



20f3RhQf7c%g,N,)N^ 
y 



'R'^if.) 



A. 



20(3RU-rn\gMN^[i2nT)' - M^,] 



(39) 



In the third and fourth lines the difference of the actions was expressed in terms of 
the radii of the circles and in terms of the glueball mass scale and the temperature, 
respectively. The difference of actions scales as A^^ in the 't Hooft large N^ limit, as 
expected. 

The difference ()39|1 between the action densities is the same, up to a factor of jS, 
as the difference between the free energy densities. It is now easy to see that for low 
temperatures, 27iR < P, the background (Hj) has a lower free energy so it dominates, 
while the opposite is the case for high temperatures, as we claimed in section 3.1. 
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